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Outline 


* Intro - codes, conventions and studies 
« Wind Convention 
* Mesh refinement 
* Accuracy requirements 

« “Axibody” - Body of revolution 


¢ “LM 1021” - Wind tunnel model of full configuration from 2014 boom workshop 


« Summary 


Introduction 


* Propagation using sSBOOM (v2.5)* for all cases 
« Augmented Burgers’ eq. 
+ Finite-difference with space-operator splitting 
* Most runs under 1 min on laptop 


* Loudness metrics computed with LCASBt 

* Applied current “best practices” 

* Mesh refinement study done for both geometries using std. atm. 
* Ran all required & optional cases 


* Rallabhandi, S. “Advanced Sonic Boom Prediction Using the Augmented Burgers Equation” 
J. of Aircraft 48:1245-1253, 2011. 


¥ Shepard & Sullivan, “Loudness Code for Asymmetric Sonic Booms( LCASB)”, NASA TP 3134, 1991 


Wind Convention 


* sBOOM uses left-handed coordinate system for wind 
* B = heading, 

* B =0° A/C pointed East 

* Clockwise = +B 


* sBOOM wind tables are in meters vs m/s 


* x and y inputs are wind components (“blows toward”) 


— = (x, y) =(1, 0) is tail wind if heading is East 
oo 1 (x, y) = (0, 1) is tail wind if heading is South 
ie 
* NX (x, y) = (1, 1) is tail wind if heading is South-East 


ty 
Net resutt is that sign on 
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Mesh Convergence 


Sensitivity of noise output to mesh refinement 


* Propagation code is solving augmented Burgers’ via finite difference 
* Need to make sure we’re getting mesh converged result 


* Mesh convergence is case dependent 
* Do for each case, assume std atmosphere 


* Dissipation due to truncation error directly impacts accuracy, resolution 
requirements are driven by need to minimize error in propagation 
* Initial signal typically has < 2 k pts 
* Propagation typically requires 20-50 kpts 


Mesh Convergence 


Sensitivity of noise output to refinement of the propagation mesh 
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Mesh Convergence 


Sensitivity of noise output to refinement of the propagation mesh 
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Mesh Convergence 


Sensitivity of noise output to refinement of the propagation mesh 


A-Weighted Loudness 


Perceived Loudness 
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+ Both dB(A) and PLGB show similar behavior 
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Time Miliseconds (ms) + Lower than 35 kHz, noise outputs drops quickly 
+ However, mesh convergence not convincing, even at higher frequencies 
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Mesh Convergence Mesh Convergence 
Sensitivity of noise output to refinement of the propagation mesh Sensitivity of noise output to refinement of the propagation mesh 
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* Oversampling introduces higher 
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* Possibility of aliasing due to 65.9 68k pts/L 500 kHz frequencies which may effect 65.9 68k pts/L 500 kHz 
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* Spline data at high resolution with Akima spline 
« Pass high-resolution data from splined signal into sSBOOM to avoid aliasing 


high-frequencies 


* Need to investigate more 
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Caveats on Accuracy Requirements 


Decibels are logarithmic units! 
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Caveats on Accuracy Requirements 


Decibels are logarithmic units! 


Ground uauianads 4 Double the loudness — ~10 dB more sensed loudness level (psycho acoustic) 
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16, 
: * Propagation accuracy requirements increase logarithmically as 
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, * Sampling frequency for a 90 dB signal is likely to be insufficient 
4 for a 70 GB signal 
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Caveats on Accuracy Requirements 


Decibels are logarithmic units! 


Ground isaaness 4 Double the loudness — ~10 dB more sensed loudness level (psycho acoustic) 
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* Propagation accuracy requirements increase logarithmically as 
signals get quieter! 


* Sampling frequency for a 90 dB signal is likely to be insufficient 
for a 70 GB signal 
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Axibody 


Shaped axisymmetric body of revolution 


Conditions: 
M.. = 1.6 
Altitude = 15849.6 m (52 kft) 
Lref = 42.98m (141 ft) 
r/L = 3.0 at signal extraction 
Ground reflection factor = 1.9 
Heading East (8 = 0°) 


Cases: 
Required: Atm #3 
Optional #1: Std. Atm. 
Optional #2: Atm #4 
Optional #4: Std. Atm. with 
70% humidity 


Axibody 
Close near field signal 
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Time Miliseconds (ms) 
* Compared 2 different closures (both linear ramps) gave consistent results 
* Closed signal using linear ramp to 0 at 100 m 


* Ground signals & noise both virtually identical 


Axibody 


Ground signature — Atmosphere #3 vs Standard Atmosphere 
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+ Atmosphere 3 was required case 
* ~0.6 PLdB louder than standard atmosphere 


Axibody 


Ground signature — Standard Atm. 


vs Atmosphere 3, ¢ = {-45°, 0°, 45°} 
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Axibody 


Ground signature — Standard Atm. 


vs Atmosphere 4, ¢ = {-45°, 0°, 45°} 
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Axibody 


Ground signature — Standard Atm. 


vs Standard Atm. + 70% Relative Humidity 


Axibody 


Raytubes, standard atmosphere 


Colored by raytube area 
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Raytubes, Atm #3 


Colored by raytube area 
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Axibody 


Raytubes, standard atmosphere 70% humidity 
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Perceived loudness at ground level 


— — IB Required [] Opt#1 [ Opt#2 [ Opt#3 
15 | | 
| | Atm #3 Std. Atm Atm #4 Std. Atm 70% RH 
| 85 
| 
> mo 80 - 
= 10 zg 
i= 
Qa oa 
® - 75 - 
$ A Z, 
3 2 7- LA 
= i= te’ 
5 ao) ey 
3 34 
oS 65 - 8! 
Fo oO gor 
5 £60 8 
moe 3 a 
= S55 - Z 
* io 
Optional #3 Crossr, 50 
ange Dis = -45° o=0° 
tance (km) x . 
Axibody Axibody 
dB(A): Optional Atmosphere profiles dB(C): Optional Atmosphere profiles 
Hi Required (| Opt#1 J Opt#2 (i Opt#3 HM Required §) Opt#1 J Opt#2 (i Opt#3 
Atm #3 Std. Atm Atm #4 Std. Atm 70% RH Atm #3 Std. Atm Atm #4 Std. Atm 70% RH 


70 


SS 


54 
fe) 
=) 
> 
o, 

gor 
51 
> 
co) 
a" 


A-Weighted Loudness, dB(A) 
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C-Weighted Loudness, dB(C) 
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Near Field Signatures 


Signals at 222.6m from flight path 
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Ground signature: Standard Atmosphere, 
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Ground signature — Atmosphere #1, , d = {-30°, 0°, 30°} 
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Ground signature — Optional #2, Atmosphere #4, , d = {-30°, 0°, 30°} 
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Ground signature — Optional #3, std/ atm. +70% relative humidity, = {-30°, 0°, 30°} 
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Standard Atmosphere, 70% Rel. Humidity 
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Slightly quieter (0.3-0.4 dB) than in std atmosphere 
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Raytubes, Standard Atmosphere 


Colored by raytube area 
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Optional #1 
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Raytubes, Atmosphere #1 
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Raytubes, Standard Atmosphere with 70% relative humidity 
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Perceived loudness at ground level 
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Signal cutoff 
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Summary 


+ Applied sBOOM and LCASSB for all required and optional cases 
+ Showed difficulty in obtaining mesh-converged loudness metrics 
+ Noted issues due to oversampling 


+ Noted that resolution requirements increase as signals get quieter 


+ Open questions & Opportunities 
+ Impact of step-size on mesh convergence? 


+ Splining of the input 


+ Consider higher-order discretization to control truncation error at 


+ Automatic output-based adaptation of sBOOM propagation 


Question’ 


understanding the ray path sorraas 


Sound 
generated 


Ray that sound 
travels along 


Effect of wind on raytube (SBOOM coords) 


here is a snapshot of the raytube with the wind. 
tube is scaled by local raytube area 
* ” 
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*“% sBOOM_2.5 -v -d” reports: 
Ground intercept for Azimuth: 0.000000 
— Propagation time: 68.010707 (s) 
— X = 51811.993216 (ft) Y= 1478.793091 (ft) 
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